Restricted quantum-classical correspondence 
and counting statistics for a coherent transition 
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The conventional probabilistic point of view implies that if a particle has a probability p to 
make a transition from one site to another site, then the average transport should be (Q) — p 
with a variance Var((5) — (1 — p)p. In the quantum mechanical context this observation becomes 
a non-trivial manifestation of restricted quantum-classical correspondence. We demonstrate this 
observation by considering the full counting statistics which is associated with a two level coherent 
transition in the context of a continuous quantum measurement process. In particular we test 
the possibility of getting a valid result for Var(Q) within the framework of the adiabatic picture, 
analyzing the simplest non-trivial example of a Landau-Zener crossing. 



I. INTRODUCTION 

The discussion of quantum-classical correspondence 
(QCC) in the mathematical physics hterature is as old 
as the history of quantum mechanics. Traditionally one 
has in mind complicated semiclassical approximation 
schemes for calculating the propagation of wavepackets. 
A more elaborated analysis of the spreading process [l| , 
leads to the distinction between robust restricted QCC 
and fragile detailed QCC. Restricted QCC means that 
(^)qm ~ (^)ci holds only for a restricted set of observ- 
ables, and survives even in the presence of diffraction, 
while detailed QCC means (^)cim ~ (^)ci for all the well 
behaved observables, and requires smooth potentials as 
in the traditional formulation. 

Controlling atoms in a few site system is state of the 
art d, Q ■ The prevailing studies are focused in restricted 
questions such as "what is the probability for a transition 
from one site to the other site" [1, Q . But what about 
the associated noise [1, 0| and the full counting statistics 
(FCS) of a quantum transition? We argue below that 
the simplest example of a single particle in a closed two 
site system can be worked out exactly so as to illuminate 
the essence of restricted QCC in the context of FCS 
studies. This has the side benefit of shedding new light 
on the vague picture in the pioneering works 8, 9] about 
shot noise, where the discussion of FCS is immersed 
in complicated diagrammatic calculations involving a 
many-body system of Fermions in an open geometry, 
hence obscuring the simple physics that underlays the 
bottom line results. While most follow-up publications 
about FCS [IS EH, El [H, m are aimed in studying 
the distribution of transmission eigenvalues, our concern 
below is with some fundamental aspects that distinguish 
quantum dynamics from its classical stochastic analog. 

We consider the simplest non-trivial example: the adi- 
abatic crossing of a particle from one site to another 
site, which involves a two level Landau-Zener transition 
[l5, 16, 17, 18]. Within the framework of a classical prob- 
abilistic point of view the particle has some probability p 
to make the transition. In any particular realization the 



particle either makes the transition or not. Accordingly 
the integrated current from the first to the second site is 
either Q = 1 or Q = respectively. It follows that 



(Q) = p 



with the variance 



Var(Q) = {l-p)p 



(1) 



(2) 



This classical probabilistic point of view does not hold 
in the quantum mechanical reality. The characterization 
of the full counting statistics requires some care in 
the description of a continuous quantum measurement. 
Accordingly one has to distinguish between the naive 
mathematical definition Q of the Q probabilit y di stribu- 
tion P(Q), and the proper physical definition [20| of the 
quasi-probability distribution 'P{Q;x), where x signifies 
the strength of the interaction with the detector. It is 
important to realize that in general P((5;a;=0) is not 
the same as P((5), but still they have the same first and 
second moments. 

This paper has two parts. In the first part we estab- 
lish restricted QCC for the counting statistics using gen- 
eral principles. In the second part we study whether the 
adiabatic approximation can be used in order to derive 
leading order results. The importance of the latter anal- 
ysis becomes apparent once one tries to obtain results for 
more complicated multiple-path geometries where exact 
solution for P((5) is out of reach. We further discuss this 
latter point in the concluding section. 



II. RESTRICTED QCC FOR MOMENTS OF Q 

For a general network that consists of several sites we 
can define a set of occupation operators A/" and current 
operators T. In the Hcisenberg picture the time deriva- 
tive of any occupation operator equals a sum over the 
ingoing current operators. In order to demonstrate the 
general idea of restricted QCC we consider below the sim- 
plest example of one particle in a two-site system. The 
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model Hamiltonian in the position basis is 
_ , at/2 c 



(3) 



where c is the hopping amphtude. We define an occupa- 
tion operator J\f and a current operator X as follows: 
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(4) 



(5) 



Using the Heisenberg picture the two are related by the 
equation 



-m = i{t) 



(6) 



where X{t) — U{tyXU{t), and U{t) is the evolution op- 
erator. Consequently it is useful to define a counting 
("transported charge") operator: 



Q= / X{t')dt' 
la 



(7) 



From the Heisenberg equation of motion Eq.([6|) it fol- 
lows that the occupation probability {Af{t)) and the 
probability current {X{t)) are related by a continuity 
equation. But in fact we can derive a stronger statement. 
By integrating Eq.([6]) over time one obtains 



Q = U{t) - Af{0) 



(8) 



Assuming that the particle is initially in the left site, 
which is a zero eigenstate of M , one concludes the non- 
trivial relation 



for fc = 1,2. 



(9) 



This relation between the fcth moment of the counting 
operator (as determined in the Heisenberg picture by 
an expectation value at the reference time t=Q) and the 
fcth moment of the occupation operator (as determined 
in the Schrodinger picture by an expectation value at 
time t) does not hold for fc > 2, because the operators in 
(0|[7V(t) - A/'(0)]''|0) are non-commuting. 

Thus for a coherent quantum transition where initially 
{M^)q — 0, we find at the end of the process {N'^)t = P, 
and consequently it follows that (Q'"') —p for k — 1,2, 
leading to Eq.([T]) and Eq.®. But the higher moments 
would be different compared with the classical expecta- 
tion, as further discussed in the next section, and there- 
fore we have here a very simple example for restricted 
rather than detailed QCC. 



III. THE FCS - NAIVE RESULTS 

Both Eq.((T]) and Eq.Q can be derived on the basis 
of a classical probabilistic point of view. The classical 
reasoning is based on the following idea regarding the 
counting statistics: 



P(Q) - { 



for Q = 
for Q = 1 



(10) 



The fcth moment of this probability distribution is 
(Q'^) =p, leading to Eq.((l|) as weU as to Eq.©. But 
we are going to explain that the classical reasoning is 
wrong: It is wrong both according to the naive math- 
ematical definition of P(Q), and also according to the 
proper physical definition of P((5; x), which we review in 
Appendix A. Consequently in the present section we de- 
rive the mathematically correct results for P{Q), while in 
the next section we work out the physically meaningful 
result for P(Q; x). 

Let us elaborate on the straightforward procedure for 
naive FCS calculation. The first step is to write the cur- 
rent operator in the Heisenberg picture: 



lit). 



{n\U{tyXU{t)\m) 



(11) 



It is important to realize that the current operator and 
hence the counting operator have a zero trace. Conse- 
quently integrating over time we get 



+Q\\ 



iQj_ 

-Q\\ 



(12) 



The first two moments (Q) and (Q^) are obtained from 
this matrix, leading to the identifications: 



(Q> = 

Var(Q) = 



\Qi 



(13) 
(14) 



The zero trace property also implies that the eigenvalues 
of Q are opposite in sign: 



Q± 



±a/(Qii)2 + IQ±P 



(15) 



On the basis of the continuity equation we can argue that 
{Q)i; l£ 1 for any preparation. Therefore we must get 
\Q±\ < 1. In fact we can deduce much more on the basis 
of restricted QCC. Observing that Eq.((l3| and Eq. fTil) 
should be in agreement with Eq.((T|) and Eq.([21) one de- 
duces that 



P 



— V {^^p)p X PhaseFactor 



(16) 
(17) 



If we could regard Q as a conventional observable, then 
upon measurement its observed values would have the 
distribution 



P(Q) 



p- 
p+ 



for Q = Q- 
for Q = 



(18) 
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This distribution is characterized by two parameters be- 
cause Q± are opposite and sign, and p± sum up to unity. 
Restricted QCC provides the two equations {Q) = p and 
(Q^) = P that can be solved, leading to 



P± 



(19) 
(20) 



where the expression for Q± is in agreement with 
Eas. (fT5)l - (|17p of the previous paragraph. A look-alike re- 
sult had been obtained for shot noise of Fermions using 
complicated diagrammatic techniques Q . An equivalent 
way to express Eq. (fT8|) is to say that the fcth moment is 



P+Q, 



P-Q- 



p 



(21) 



where fc = 0,l,2,3,... and [...J stands for the integer part 
(i.e. rounded downwards). The corresponding classical 
result is p irrespective of k. One may say that we have 
encountered here the simplest example of restricted QCC. 



IV. THE FCS - PHYSICAL RESULTS 

In a later publication Q the naive mathematical def- 
inition of the full counting statistics has been criticized. 
The most illuminating approach [20| is to analyze the re- 
duced dynamics of the detector. Using the standard von- 
Neumann pointer scheme [2l| and transforming to the 
Wigner representation one deduces that the final state 
p{q, x) of the pointer is the convolution of its initial state 
with a kernel P{q—q'; x), as defined in Eq. (jA5|) . It follows 
that 



P(g;a;=0) = 
1 

2^ 



(22) 



7-e-«('-/2)Q 



-"^''dr 



can be regarded as a quasi-distribution that describes 
the full quantum statistics. It is of course a physically 
measurable object. The naive mathematical definition is 
obtained if we ignore time ordering: 



P(0) = TT- 



1 
271 

m Q)) 



(23) 



In general the calculation of P(Q; a;=0) is very compli- 
cated. But we can gain some insight by considering the 
simplest case of a Bloch transition, which is Eq.© with 
a — 0. After time t the probability to find the particle 
in the second site is p = [sin(c<:)]^. The expectation value 
in the integrand of Eq. is 



P(r) = 1 - 



(r/2) 



sm 



ct^l + (r/2)' 



(24) 



[More conveniently one can use instead of Eg. ip^ the 
equivalent expression Eq. (jA5|) . where the Us can be in- 
terpreted as spin 1/2 rotation matrices]. In order to get 



P((5; a;=0) we have to Fourier transform (FT) the func- 
tion P(r). For the purpose of discussion let us assume 
long times {ct » 1) so as to have separation of scales. 
Then we can distinguish between a central part where 
P(r) « 1 -|- ipr — (l/2)pr^ -|- ... as implies by restricted 
QCC, and oscillatory far tails where P(r) « cos(ctr). The 
FT of the central part gives a non-singular exponential- 
like piece, while the FT of the tails contributes two delta 
functions {1/2)S{Q ± ct) which are screened by negative 
clouds as illustrated in Fig. 1. We can regard the singu- 
lar delta functions at Q = zLct as the remnants of those 
that are centered at Q± = ^^/p in the naive calculation. 
It is interesting that for long time their location reflects 
the eigenvalues of the current operator rather than the 
eigenvalues of the counting operator. 

Fig. 2 is a caricature that illustrates the significance of 
the FCS results with regard to the outcome of a quantum 
measurement. We sketch the final probability distribu- 
tion of the Von Neumann pointer as implied by different 
"conceptions" . The upper panel is based on the classi- 
cal expectation of having either Q = or Q = 1 displace- 
ment. The middle panel is the naive quantum mechani- 
cal prediction that suggests Q = Q± displacements, while 
the lower panel is the actual quantum prediction based 
on a convolution with P((5; x). An actual numerical illus- 
tration for the outcome of this convolution is presented 
in Fig. 3. Needless to say that in principle one can mea- 
sure not only the q distribution but the whole probability 
matrix so as to determine P((5; x) and in particular the 
quasi probability distribution P{Q; x=0) via a deconvolu- 
tion procedure. But looking in Fig. 3, it is nice to realize 
that the main features of the FCS are not smeared and 
pop to the eyes even without any deconvolution. 



V. FCS IN THE ADIABATIC 
APPROXIMATION 

In the following sections we would like to examine the 
capabilities of the leading order adiabatic approximation 
in obtaining physically significant results for Var(Q). We 
consider the model system of Eq.Q. Initially (at t= — oo) 
the particle is in the ground state, which is the left site. 
The probability at t=oo to find the particle in the right 
site is 



where 



p 


= 1- 
















\ o C^l 




TT " 


exp 


-2tt— 


= exp 






a 







(25) 



(26) 



is the Landau-Zener probability to make a transition 
from the lower to the upper energy level. The rate of 
the driving is characterized by the dimensionless param- 
eter 
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(2c)Va. 



(27) 
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The straightforward way to calculate Var(Q) is also the 
most complicated one, because it requires an explicit 
evaluation of the evolution operator U{t). For our two 
site system this can be done using parabolic cylinder 
functions p^ . But having in mind more complex sys- 
tems [2^, for which exact solution arc not available, we 
would like to make the explicit calculation within the 
framework of the adiabatic approximation: 



u{t) « 51 hw) / 



En{t')dt' 



n{to) 



In this expression \n{t)) and En{t) are the so-called adi- 
abatic states and the adiabatic energies The as- 
sociated expression for the matrix elements of the time 
dependent current operator in the Heisenberg picture is: 



n m) 
n{t)\I\m{t)) X exp 
ice**(* 



to 



where 



m = 



yJiatf + (2c)2 dt' 



(28) 



(29) 



Using the zero order {0{a^)) adiabatic states we get for 
the off diagonal terms ~ Qlz where 



cosh(z)e**(^)dz 



(30) 



Above we use the substitution t = {2c/a)T followed by 
r = sinh(z). 

Before we go on with the calculation of Qlz we would 
like to comment on the calculation of the diagonal terms 
in Eq. (P5)) . This terms are trivially related via Eq. (fTB)) to 
the Landau-Zener transition probability Plz of Eg. ipS)) . 
The leading order estimate of this probability is based on 
a conventional time dependent treatment 17], leading to 



1 

2„_ 



1 



2 cosh(z) 



(31) 



It is important to realize that in the zero order adia- 
batic approximation we get zero for the diagonal terms 
of Eq. ([28|) . because the zero order adiabatic states are 
time-reversal symmetric and hence do not support non- 
zero average current. If we use the first order {0{a^)) 
adiabatic states we get for the diagonal terms ±1 but 
still miss the non-adiabatic Plz correction. It is there- 
fore non-trivial and requires verification that Q± w Qlz 
is a valid approximation. 



VI. THE CALCULATION OF Qlz 



In order to evaluate Qlz of Eq. ([30)1 we use the contour 
integration method as in Ref. [l3| ■ The explicit expression 
for the phase $ as a function oi z = x + iy is 



+ ^ sinh(2z) 



(32) 



X -f - sinh(2x) cos(2y) 



+i \ y + - cosh(2a;) sin(2y) 



The contour of integration in Eq. (|30p is y — 0, but we 
would like to deform it into the complex plane so as to get 
rid of the rapid oscillations of the phase factor, and have 
instead a smooth monotonic variation. The deformed 
contour is displayed in Fig. 4. The phase is pure imagi- 
nary along the curves C_ and C+. At zq ~ + i(7r/2) we 
have $ = «(7r/4)7, while coshz w i(z — zo). Consequently 
in the Plz integral of Eq. (PT|) we have a pole leading to 
the standard LZ results (disregarding the prefactor which 
requires higher orders). But in the case of the Qlz in- 
tegral Eg. ((30)) we do not have a pole: rather we have 
to consider the non singular part that comes from the 
integration along the C± curves. One observes that 



.dz = 2n ...dz 

C-+C+ JC+ 



fix) e**(^)da; 



where 



i^{x) 



TT — arccos 



2x 



1 - 



2x 



sinh(2a;) 
1/2 



sinh(2x) 



cosh(2a;) 



and 
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V2 



2x 






-1/2 ^ 


sinh(2a;) 





cosh(4a;) — 1 



sinh(a;) 



Deep in the adiabatic regime (7 3> 1) the integration is 
dominated by the small x interval where f{x) cx 7a; and 
[$(a;)— $(0)] cx jx^. Accordingly 



Ql 



7^/^ exp 



(33) 



Inspired by the calculation procedure of Plz in Ref.[l6|, 
our speculation is that also here the pre-exponential term 
would be renormalized to unity by higher orders, which 
would imply consistency with the expected (exact) result 
|Q±| = V(l - ^lz)Plz that follows from Eq.fTTl). 
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VII. DISCUSSION 

We have demonstrated how restricted QCC can be es- 
tabhshed and utihzed in order to determine the counting 
statistics of a quantum coherent transition. It is impor- 
tant to reahze that this procedure has some hmitations. 
Namely, if we had considered not a two site system but a 
more complex system with multiple path geometry, then 
the same considerations would not allow to deduce the 
variance Var(Q) of the integrated current. In particu- 
lar one wonders what happens if a particle has two op- 
tion al p aths available to get from one site to another 
site ^2]. Within the framework of the classical proba- 
bilistic theory, a splitter would imply a noisy outgoing 
current. But in Ref.[2^ we argue that the coherent split- 
ting of a wavepaket is not noisy, and furthermore that 
the whole study of fluctuations in quantum stirring de- 
vices requires to go beyond QCC considerations. Con- 
sequently the problem that has been raised in Ref. j23| 
has motivated us to study the capabilities of the leading 
order adiabatic approximation for the purpose of calcu- 
lating Var(Q). 

The calculation of the Landau-Zener transition proba- 
bility Plz and the associated dispersion Qlz in the leading 
order adiabatic approximation yields a contour integral 
in the complex plane. While the integral is dom- 
inated by a pole, the Qlz integration is related to the 
corresponding principal part. Restricted QCC implies 
that Plz and Qlz are related. Our analysis has demon- 
strated that the dominant exponential term is correctly 
reproduced, but not the pre-exponential term which ap- 
parently requires infinite order. 

The traditional QCC principle is based on the semi- 
classical approximation and implies detailed QCC for all 
the moments up to the quantum resolution limit. In 
contrast to that restricted QCC is very robust, and sur- 
vives even in the presence of diffraction [l| . The Landau- 
Zener problem is possibly the simplest non-trivial exam- 
ple where this idea can be demonstrated. For some more 
complicated (chaotic) systems restricted QCC can be es- 
tablished as an approximation in the long time limit, on 
the basis of the short time perturbative analysis which 
is extrapolated using the central limit theorem. It would 
be interesting to explore the implications of coupling to 
the environment in this context 1241 . 



Acknovifledgments — The original motivation to deal 
with counting statistics in the context of closed systems 
came from a discussion with Oded Agam (HUJI). Yuly 
Nazarov (Delft) is acknowledged for highlighting some 
major observations regarding full counting statistics fol- 
lowing Refs.[8l.[9l[20j. DC thanks Miriam Blaauboer for a 
most fruitful visit in Delft TU, during which this work has 
been crystallized. We also enjoyed talking with Yaroslav 
Blanter (Delft), Yuval Gefen (Weizmann), and Shmuel 
Fishman (Technion). The research was supported by 



grants from the Deutsch-Israelische Projektkooperation 
(DIP), and the US-Israel Binational Science Foundation 
(BSF). 



APPENDIX A: THE VON-NEUMANN 
MEASUREMENT SCHEME 

In this section we present a short simple derivation of 
the main result of Ref.jl^] regarding the measurement 
of the full quantum statistics. The original derivation 
has been based on an over-complicated path integral ap- 
proach. 

The coupling of the system to a Von-Nuemann pointer 
PH whose canonical coordinates are (i, q) is described 
by the Hamiltonian 



Wtotai = n{t)~ IX 



(Al) 



The states of system can be expanded in some basis \n), 
and accordingly for the system with the detector we can 
use the basis \n,x). The representation of the evolution 
operator is 



U{n,x\no,xo) 



(A2) 



where U[x] is a system operator that depends on the 
constant parameter x. We formally write its explicit ex- 
pression both in the Schrodinger picture and also in the 
interaction picture using time ordered exponentiation: 



U[x] = T exp 



i [ {n-xl)dt' 
Jo 



U[0]Texp 



IX I I{t')dt' 





(A3) 



The time evolution of the detector is described by its 
reduced probability matrix 



pt{x",x') 



E 

U{n, x'^)Uin, x'\n'^, x',)* fQ.^Poijr'i. ^^o) 



Pq{x",x') 



where p"^'' = |V')(?A| is the initial state of the system, 
and pq(x'^,x'q) is the initial preparation of the detector. 
Transforming to the Wigner function representation we 
get the convolution 



Pt{q-,x)= P{q~ q';x)po{q',x)dq' 



(A4) 



where 



P{Q;x) 
1 

2^ 



(A5) 



U[x - (r/2)]t[/[a; + (r/2)] t/;^ e-''^''dr 
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Fig.l: Numerical calculation of the full counting statistics kernel 
P(Q; x=0) for a coherent Bloch transition from one site to another 
site. The evolution is governed by the Hamiltonian of Eq. l[3j with 
= and ct = 5. The thick bars represent delta functions. See 
the text for further explanations. 



Fig. 3: Numerical calculation that illustrates the final probabil- 
ity distribution of the Von Neumann pointer q corresponding to 
the lower panel in Fig. 2. The parameters are the same as in 
Fig. 1. The initial preparation is a Gaussian wavepacket centered 
at 5 = a; = of width = 1.2 (solid line), and = 0.8 (dashed 
line), and = 0.5 (dotted line). Its evolution has been calculated 
using Ea.l lA4l l with Eg. llASI l. 
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Fig. 2: A caricature that illustrates the final probability distribu- 
tion of the Von Neumann pointer q. It is assumed that the initial 
state is a wavepacked that is concentrated at g ~ 0. The upper 
panel is based on the classical expectation of having either Q = 
or Q = 1 displacement. The middle panel is the naive quantum 
mechanical prediction that suggests Q = Q± displacements, while 
the lower panel is the actual quantum prediction based on a con- 
volution with P{Q;x). 
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Fig. 4: The calculation of Plz and Qlz is done by deforming the in- 
tegration contour into the complex plane. The integration is along 
the curves C± where the integrand is non-oscillatory, and along the 
arc C'o that encircles the pole P at zg = + i{iT/2). The contri- 
butions of the (non-displayed) vertical segments that connect C— 
and C+ to the real axis at infinity can be neglected. 



